ITEP-TH-64-99 



Dirac operator spectral density 
and low energy sum rules 



K.ZyablyukQ 

O ' Institute of Theoretical and Experimental Physics 

o ■ 
o 

(N 



^ ■ B.Cheremushkinskaya 25, Moscow 117259, Russia 



Abstract 

The spectral density of euclidean Dirac operator is investigated in partially quenched 
CN ', QCD with arbitrary quark masses. A representation of scalar and pseudoscalar cor- 

relators in terms of the spectral density is discussed. The spectral density, obtained 
^T) • from the partially quenched chiral perturbation theory is shown to be compatible with 

K*" ! low energy sum rules, obtained earlier. 

o 
o 

1 Introduction 

^ ■ One of most important characteristics in QCD is the spectral density of Dirac operator 
^ 1 

9"'- PW = T7<T.^i^'^n[A])>A, (1) 

Oh! ^ 

' where A„[yl] are eingenvalues of the Dirac operator D with gluon field A in 4-dimensional 
euclidean volume V: iDun = XnUn- The average < . . . >a denotes integration over gluon 
^ ! fields with Yang-Mills action and the product of fermionic determinants det {D + mj). 
^ I Due to this averaging the spectral density depends on the quark masses mj. 

A useful tool for the investigation of the spectral density is partially quenched QCD 0. 
In addition to usual quark with masses rrii it includes one unphysical (valence) quark with 
mass m^i and one unphysical quark of opposite statistics with mass m„2- The generation 
functional of this SU{N + 1|1) theory is written as follows: 

N 

I I H-"! , I / / —I— III-, I 

ZP'^{mi, thn; m^i, m^a) 



Then one finds the condensate of unphysical fermionic quark for equal masses m^i and my2'- 



(3) 
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It depends on the masses of all quarks, but we shall write explicitly only the dependence on 
the unphysical quark mass If we put it equal to the mass of some physical quark, we 
get its chiral condensate: 

S(mi) = - < QiQi > (4) 
The spectral density is the discontinuity of S(m„) across the imaginary axis: 

27rp(A) = S(iA + 0) - S(iA - 0) , (5) 

This equation with condensate E(m) given by (0, ^ is just trivial rewritting of the definition 
(0) of the spectral density. Nevertheless it allows us to compute the spectral density using 
well known rules and methods of evaluating the Feynman graphs. 

Indeed, the condensate of unphysical fermionic quark E(m„) is the fermion Green function 
at equal coordinates. It can be obtained from the quark propagator S{p) in momentum space 
in d dimensions as the expansion over QCD coupling constant a^: 



(27r)' 

d'^p 12m , f d'^p d'^k 2id - 2)ipk) + id - A)p^ + dm'^ 
647rn' m i — — 1- Oin ) 

^ ,Vil-d/2) 7 1) r2(2 - rf/2) 

(4^)d/2 ^vra.m _ ^^^^^ + U[aj , [b) 

"tr" here is the trace over color and spinor indices, tr 1 = 3 ■ 4; euclidean space is assumed. 
The two-loop integral in can be computed with help of technique described, for instance, 
in . The result (|D is divergent in four dimensions. If we put o? = 4 — 2e and expand over 
e, we get the following result: 



f 3 , cts 

L m) = — < const + -— In — + — 



--2Vn^-3 1n< 



(7) 



where /x^ = Ane^'^^^'^^fi'^ according to MS convention. The terms ~ are not important 
for our purposes, since they are regular for < 0, while logarithmic terms give us desired 
discontinuity. The term with 1/e is removed by renormalization of the mass and condensate 
itself. Now we may apply the equation (|) in order to get the spectral density with as 
correction: 

3IA 



P(A) 



47r2 



1 - 




(8) 



Notice, that the spectral density in QCD depends on the quark masses only starting from 
level, where the loops with usual quarks appear. 

The equation @ works in the QCD range of validity, i.e. for |A| ^ Aqcd- In the region 
|A| ~ Aqcd the formula is useless, since nonperturbative effects become important. Some 
attempts to extract the behavior of the spectral density in the region |A| ^ Aqcd from low 
energy sum rules were undertaken in [Q, Two-point scalar correlators, known from 

the chiral perturbation theory (ChPT), can be presented as certain integrals over the spectral 
density. This analysis demonstrates linear behavior p(A) ~ {N"^ — 4)|A| for |A| 3> m. 
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The spectral density has been obtained in ||^ within 1-loop partially quenched chiral 
perturbation theory |^ with quarks of equal mass. This result correctly reproduces the 
3-point scalar correlator, obtained in Q]. 

In section 2 of this paper the result of is extended to the case of arbitrary quark 
masses. We take into account complete 1-loop generation functional of partially quenched 
ChPT with 0{p^) terms, which absorb divergencies. In section 3 the scalar and pseudoscalar 
correlators are analyzed with help of spectral density representation. Many of them can be 
reduced to the condensate S(m^,) (and its derivatives) and to the correlators of topological 
charge densities. Our result is shown to be in agreement with the sum rules, obtained in P, 
i- 



2 Spectral density in partially quenched ChPT 

At low energies the chiral perturbation theory [||, ^ is a good replacement of QCD. Its 
degrees of freedom are N"^ — 1 mesons instead of quarks. In analogy with this construc- 
tion one supposes that the low energy approximation of the partially quenched QCD (Q) is 
partially quenched ChPT with chiral symmetry SU{N + 1|1) 0. The effective degrees of 
freedom are mesons (pij of the same statistics as the product of relevant quarks qigj. 

The generalization of the SU{N) chiral lagrangian to the partially quenched case is 
obvious: everywhere the flavor trace should be replaced with supertrace. Nevertheless, it is 
simpler to consider the U{N + theory with super-?]' instead of pure SU{N + 1|1). In this 
case there is no need to impose constraint Str0 = and one can consider diagonal fields 0ii 
as independent variables. In the final results one can take the limit of heavy super-?]' which 
corresponds to the SU{N + 1|1) case. 

The following terms of the chiral lagrangian are responsible for the chiral condensate: 

L = -K)(0o) + V,{<l)o)StT{d,U^d^U) + 2ReStr [VMmU] + Vsl^o) d,<Podyo (9) 
where U is unitary matrix of meson fields: 

U = exp y-^r- j 5 00 = Str(/) = 0n , = 

m = diag(mi, . . . , m^r, m^i; mi,2) is the quark mass matrix, F is pion decay constant. The 
freedom of field redefinition (pij (pij + Sijf{(j)o) can be used either to vanish V3 or to make 
V2 real. We shall use the last choice, since in this case the calculation is simpler. Expanding 
the lagrangian (|D up to quadratic terms, one needs the following constants: 

vm = -J , v,{o) = , K,"(o) = /i^ , ^3(0) = - . 

There is also term V2"(0) Strm0Q. It belongs to the same class, as the terms of order O(p^), 
so the constant V^"(0) is supposed to be small. Since it contributes to the condensate only 
at the loop level, we shall ignore it. 




N + 1 



N + 2 
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The derivation of the one-loop chiral condensate (j^) is standard. The propagator of the 
diagonal mesons (pa was found in [Q: 



G 



+ E 



-1 



„2 _ A/f 2 



(10) 



where = 2Bmi. Later we consider only the limit —>■ oo. In this case the super-?]' 
decouples leaving us with pure SU{N + 1|1) theory. In euclidean space p"^ changes the sign 
and the chiral condensate (E) can be written as 



N 

E 



1 



1 



N 



-1' 



~{ p'^ + B{mi + m^) (p2 + 2Bm^Y p'^ + 2Bmk J 



When the volume goes to infinity the sum over momenta in (|TTp should be replaced with the 
integral over d'^p / {2Tr)^ . The integral is quadratically divergent, so we shall use dimensional 
regularization. We separate divergent terms according to usual conventions of ChPT 



N 



S(m^) = S^(m^) - 2B^c ^ 

1 



where 



167r' 



i=l 
d-A 



(12) 



L,-^4(r'(l) + ln4. + l) 



E''(m^) is regular in c? = 4 part. The divergencies are absorbed by the coupling constants of 
the O(p^) lagrangian. The terms without derivatives of this lagrangian: 



L(^) = Le Str (X+f/ + U+X)^ + Lj Str (X+f/ - U+Xf 
+ Lg Str (X+t/X+t/ + f/+Xf/+X) + H2 Str (X+X) 



(13) 



where X = 2Bm. Adding the contribution of these terms to (0), we get the final result for 
the one-loop chiral condensate: 



BF'- 

N 

-E 

i=l 

+ %B' 



B' 



167r2 



X dx 



{x + m^y 



f N 

E 



1 



TV 



X + TTli 



(rrii + m^) In ■ 



Birrii + m„ 



-^5](x + m,) 

1=1 

2Bm,, 



+ jV^K - 2m,) In — ^ + — (m, 



N 



i=l 



(14) 



where L'^ifi), H^{fi) are renormalized constants which depend on the normalization scale fi 



N'^ + 2 _ 4 - A 

LW=L,~^^c, Llif^) = Ls - , h;{i^) = H, - . (15) 



16A^ 



In the integral we subtracted two terms, divergent at x = p^ /2B — > cxd, so the total integral 
is finite. Due to boson-fermion loop cancellation the couplings of SU{N -l- 1|1) partially 
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quenched ChPT are renormalized exactly in the same way as the ones of usual SU{N) 
theory. Obviously, the total result (llip does not depend on /i. 

Now one can compute the discontinuity of the condensate S(m^) across the imaginary 
axis according to (1) to find the spectral density. Nevertheless, the r.h.s. of the equation 
(0), as it stands, does not have appropriate analytical properties for complex rriv. Indeed, 
according to Banks-Casher relation the condensate must have the form E(m) = m/(m^), 
where the function f{rn?) has a cut in m? plane along half-axis m? < 0. The r.h.s. of ([T^) 
does not satisfy this requirement. Consequently, in order to get the spectral density from 
(0), the formula should be apphed with the following modification [0: 



2np{\) = Disc 



my 



S(iA) + S(-iA) 



(16) 



Notice, that in QCD the equations (|) and ([TBD give the same result (||). 
Final result for the Dirac operator spectral density: 



P(A) 



TT 



N 



-E 

i=l 



327r3 



X[X 



+ A2)2 



^ N 

E 



-1 



N 



— ^(x + m, 

i=i 



dx 



, B\X^ + m^) A 4 , 2B\X\ An 
rrii in — zAarctg — + -^^^rriim h ^tt^t |A| 



Ar2 



+ jV^"^^ ~ (327r)2L;:(/i)mi 



;i7) 



In particular case of equal quark masses mj, the integral in ([TTD disappears and our result 
coincides with eq. (83) of ref. |J7[ up to the terms, linear in mass. In our case there is no 
ambiguity among these terms, since we have included all the 0{p^) contributions. 

In deriving the eq. (0) we neglected the contribution of higher loops. Each diagram in 
ChPT brings extra factor together with loop factor {Att)^"^. Consequently the expression 
(0) for the spectral density is valid for |A|, rrii ^ AnF, which is typical expansion parameter 
in ChPT. 

Two flavor case deserves special attention. For N = 2 there appears extra term of order 
O(p^) in the chiral lagrangian: 



h 
16 



Str (X+f/ + U+xy 



16 



Str {X^U-U+xf 



+ ^{hi + h) Str (X+X) + i (/ii - hs) Re Sdet X 



the couplings k, hi are the constants of the two flavor chiral lagrangian M. In this case the 



last line of eq. (|T4D should be replaced with: 



... + B' 



213(12) (mi + 1712) + 2{hi + h3)my + [hi - 



mim2 



rrir 



The spectral density for = 2 is given by ([TtD with obvious replacement Lq /3/I6. Again, 
the high energy constants hi do not contribute to the spectral density. 
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Few steps further can be done to evaluate the integral in (0). One may notice, that the 
expression inside this integral can be written in the following form: 



f N 

E 

^1=1 



-1 



X + rrii 



1 ^ a- 



(19) 



The masses rhj are the roots (with opposite sign) of the polinomial, which stands in denom- 
inator of the l.h.s. in equation ([T9|). The constants aj have dimension (mass)^; they vanish 
when all rrii are equal. In particular case of two quarks: 



mi 



(mi + m2 



ai 



mi - 1712 



For = 3 the constants mj and aj also can be found in explicit form, but we shall not write 
them here. Knowing these constants, the integral in ([17D can be easily computed: 



/t^^tpf iS^+^j "^S^""^ 



rrii 



dx 



N~-l 



E 



ajnij 



{ (A2 + m2)2 _ 



(A — m.) In — A — m- + vrlAlm,- 
■'mi 



(20) 



3 Low energy sum rules 

The low energy sum rules are formulated for the correlators of charge densities which can be 
written as certain integrals over the spectral density. In particular, the two point correlators 
of scalar and pseudoscalar charge densities 



were expressed in terms of the spectral density in 0. This technique can be extended to the 
case of different quark masses: 



dx < P''{x)P^\<d)^ > = 25'^5^ 



J (A2 + m2)(A2 + m2) 

(A^ + minij) p(A) 
(A2 + m2)(A2 + m2) 



d\ , k^l; 
dX + 5'^''^ <u^> 



(21) 
(22) 



where < z/^ > is the topological susceptibility: 

i dx < uj{x) cu(0) > , 



■tre(G^.G'^'^) 



(23) 



Correlators of diagonal scalar charges include additional contributions of disconnected graphs. 
The insertion of the operator —i J dxS^^{x) corresponds to the differentiation over the mass 
mj. In particular 



i dx< S''{x)S^^{(}y > 



drrii 



(24) 



Since the condensate is itself the derivative, Lh.s. is symmetric over ij. 

There is no need to compute the integrals in (2TJ2^), since they can be expressed in terms 
of the chiral condensate: 

i[dx< S'^{x)S''\0y > = SK) - S(m,) ^ ^ _^ ^. ^ _^ ^. ^25) 

J TTii — rrij 

ijdxK P'^{x)P'\oy > = 6V + ^K) + s'js'^i < ^2 ^ (26) 

These equations directly follow from the limit g'^ — in the appropriate vector/axial corre- 
lators q'^q'^ J e^'^^dx < V^{x)V^{0) >. Nevertheless, the equation (pSj) within the framework 
of the partially quenched theory can be also apphed if some quarks have equal masses: the 
ratio should be replaced with derivative. The result of this limit rrii rrij is however not 
the same as the l.h.s. of ( p4|) due to different order of operations: one should compute the 
derivative dT,{my) / drrit; at first and substitute = rrii after then. 

In the same way other n-point off-diagonal correlators can be expressed in terms of the 
chiral condensate. In particular 



,-n-l 



dXi... dXn-1 < S'^^'iXi) . . . ^^"-^^"-l(x„-l)S^"^"(0) > 



[6''^^ . . . + perm{j)) ^ — , ^ for any c . (27) 

a=i n {m, - mi J 



*6 

6=1, 



where perm{j) denotes all permutations of indices j. 

A large set of sum rules can be formulated using spectral density representation. Below 
we consider the simplest realistic model with N = 2 flavors of equal mass m„ = = m. 
The charges are combined in triplet S"" = qa°'q and singlet = qq. The following sum rules 
have been formulated earlier 0,0: 



ijdx< S''{x)S\0) - 6^^P%x)P\0) > = 25"^ 



drriy 



S(m) 



2 < > 



m 



ijdx< r^5°(x)5°(0) - P^{x)P\0) > = 25"'' 



6''^8BHj ■ {2t 
dTj{m) S(m)\ 



dm m 



. f_!^ _ m , (29) 

We have expressed the integrals over the spectral density in terms of the condensate according 
to ( pi| - PBD . Equation ( PD| ) is satisfied for S(m^,) given by (0). The sum rule ( P^ gives the 
following result for the topological susceptibility: 

<v'>= -BF'm + BV[2ll{^^) -2h - (30) 

This result can be obtained directly from the two flavor chiral lagrangian [Q, which verifies 
the partially quenched approach used here. 
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